Preprint typeset in JHEP style - HYPER VERSION 



Three Roads to Probe-Brane Superconductivity 



i Bum-Hoon Lee and Matthias C. Wapler 

Center for Quantum Spacetime, Sogang University, Seoul, Korea 

(N 

i-O E-mail: bhl@sogeing.ac.kr ; wapler@sogang.ac.kr 

D 

00 

Abstract: We study a defect system of two parallel D5 probe branes in a large- N c D3 
background. Using the non-abelian DBI action, we study three different fields that can 
give rise to a superconducting phase transition: A vector (p-wave), a scalar corresponding 
to a non-trivial "separation" of the branes in the (3+1) field theory directions and a scalar 
(]} corresponding to a separation in the "internal" S" 5 (both s-wave). 

i Comparing these phases first in the a' 2 expansion, we find that the internal scalar has the 

largest critical temperature and is always thermodynamically preferred. Further, there is 
an interesting attractor behavior. 

Taking a simplified version of the full DBI action that preserves its regularity and geometry, 
we find that the divergences of the a' 2 expansion are resolved and some second order tran- 
'■O sitions turn into first order ones. In addition to some other changes of the phase diagram 

due to the structure of the DBI action, we observe that the ground state degeneracy of the 
unbroken theory is lifted. We also isolate the unphysical artifacts of our simplification. 
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1. Introduction 

In recent years, the ADS/CFT correspondence [1, 2, 3] has become a powerful tool to 
study various properties of the strong coupling limit of conformal field theories, with ap- 
plications to QCD and also more recently to some aspects of condensed matter physics. 
In principle, one has to differentiate between top-down setups that are constructed within 
string theory and imply consistency - and bottom-up setups, where the gravitational duals 
are constructed from a phenomenological point of view. In this paper, we use the former 
approach as we would like to explore what happens to a particular consistent theory. 
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There are two ways to motivate the present work: In condensed matter applications, 
there has been particular interest and a significant amount of activity in 2+1 dimensional 
systems that can be typically constructed using M2 branes [4], or as a defect in a 3+1 
dimensional background using D3-D5 [5, 6, 7, 8, 9, 10] and D3-D7 [7, 8, 11, 12, 13, 14] 
brane intersections - and also as bottom- up setups in various contexts. As our world is 3+1 
dimensional, the defect setup may be more realistic, even though there are some difficulties 
involving the stability of the D3-D7 systems [8, 7, 12]. These are based on a common way 
how to introduce fundamental matter in ADS/CFT: probe brane configurations, where one 
considers a planar AdS blackhole background. For example one considers the well-known 
AdSz, x S 5 (planar black hole) solution (above the deconfinement phase transition) from the 
decoupling limit of a stack of N c S> 1 D3 branes, that is dual to a (thermal) SU (N c ) J\f = 4 
supersymmetric Yang-Mills theory [2]. Then one inserts Nf <C N c intersecting Dp branes, 
giving Nf families of (charged) fields in the fundamental representation of the SU(N C ), 
living along the directions of the intersection [15, 16] - such as the above-mentioned D3- 
Dp intersections. 

One interesting question to ask then is what happens if one takes two (or more) D5 
probe branes of the (2+1) defect system and studies new effects that appear if they are 
physically separated in the third (3+1) direction. This is the question that evolved into 
the present work. 

Another motivation are the constructions of gravitational configurations whose duals 
show spontaneous second-order symmetry breaking phase transitions and other properties 
reminiscent of superconductors [17, 18, 19, 20]. This sparked great interest and there has 
been a significant amount of activity in this field studying effects such as the Meissner 
effect and magnetic vortices (e.g. [21, 22, 23]), Fermi effects (e.g. [24, 25, 26]) or higher 
order /backreaction effects (e.g. [27, 28]). 

Most of these constructions were however bottom-up models. While there has been 
also some activity on constructing duals of superconductors in supergravity theories that 
can be obtained from M theory [29, 30, 31, 32] with many interesting results, e.g. the 
"landscape of superconductors" [33]. In contrast, there has been only very limited activity 
- which gave rise so far only to p-wave superconductors [34, 35, 36, 26] - on constructing 
D-brane configurations that give rise to superconductivity more directly from string theory, 
give very precise information on the details of the field theory dual and give always regular 
solutions. 

In this paper, we consider two parallel defect probe D5 branes in a D3 background. In 
the field theory side, they correspond to a 2+1 dimensional system of two species of parti- 
cles, e.g. a multi-layered system as in a cuprate superconductor or multi-layered graphene. 
One way to get some interesting physics from their separation in the flat directions is then 
to turn on a non-trivial SU(2) potential in the world-volume theory. This corresponds e.g. 
to an isospin chemical potential/density [37] or any other non-trivial combination of the 
densities of the different species. Then, we can study the phase diagram and look e.g. for 
a spontaneous symmetry breaking phase transition. In addition to the separation in the 
(3+1) physical directions, we will also consider a separation in the "internal" space and an 
SU(2) magnetic field, such that we can compare the different phases and have for the first 
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time a system in which different s-wave and p-wave superconducting condensates compete 
with each other. 

We will consider both the quadratic expansion of the world-volume action and an 
adaption of the full action that allows us to study a slightly modified version of the full 
non-perturbative action - giving us e.g. consistent regular physics in the low-temperature 
limit. We will see how the non-perturbative effects change the phase diagram, and isolate 
and discuss an artifact that comes from our simplification. 

The paper is organized as follows: In section 2, we review and describe the details 
of the string theory configuration with its field theory dual. In section 3, we obtain and 
discuss the results of the different phases in the perturbative treatment of the DBI action. 
In section 4, we describe a small modification to the action that allows us to study non- 
perturbative effects, apply this to the different phases and discuss the changes compared 
to the perturbative case. Finally in section 5, we conclude and discuss the results and give 
some outlook for possible future extensions of this work. For completeness, we present 
some results related to the grand canonical ensemble in appendix A. 



2. Setup 

We start with the supergravity background of a planar black hole in AdS$, 

^ = J ("(I " rt/r^ + d4) + £ (Y^p + , < z = (2.1) 

This corresponds to the decoupling limit of N c black D3-branes dual to N = 4 SU(N C ) 
super- Yang-Mills theory at finite temperature T, living along the flat directions of the AdS 
spacetime [3, 2]. The temperature is given by the Hawking temperature T = and the 
Yang-Mills coupling is related to the string coupling g s by g^ m = ^g s . Since the AdS 
length L is given in terms of the string coupling and string length l s as L 4 = 4tt g s N c l 4 , 
the 't Hooft coupling A = Qy m N c can be written as A = -pr. Hence the "supergravity 
limit" L > l s in which the type IIB supergravity action and the solution (2.1) are valid 
corresponds to strong coupling A 3> 1. 

In practice, however, we will use coordinates that are made dimensionless with factors 
of y- = Ap, denoted by (•), such as t := ^4 = t i^T), and the inverse dimensionless radius 
u : = giving us 

T 2 / rlt, 2 N 



u 2 



In this setup all the fields transform in the adjoint representation of the SU(N C ). 
However in QCD or condensed matter physics, one also needs to consider matter that is 
charged under this symmetry, i.e. that transforms in the fundamental representation. To 
introduce the fundamental matter one then creates an intersection of "probe" Dp branes 
with the D3 branes, such that in the string theory side there are massless fields in the 
effective field theory at the intersection. From the point of view of the probe branes, 
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they correspond to endpoints of D3-Dp strings, and in the gravity side they correspond to 
fundamental fields in the (defect) field theory [38, 39]. 

Here, we use the well- known D3-D5 defect setup (see e.g. [5, 6, 40]): 



12 3 

t x y z 

background : D3 x x x x 
probe : D5 x x x 



5 6 7 8 9 



x x 



(2.3) 



The dual field theory is now the SYM gauge theory coupled to Nf fundamental hyper- 
multiplets, which are confined to a (2+l)-dimensional defect. This construction is still 
supersymmetric, but the supersymmetry has been reduced from N = 4 to N = 2 by 
the introduction of the defect. In the limit Nf <C N c , the D5-branes may be treated as 
probes in the supergravity background, i.e. we may ignore their gravitational back-reaction 
[38, 39]. 

In contrast to [8, 7, 41, 42], we want to use the non-abelian structure of the U(Nf) 
theory for the simplest case of Nf = 2. Furthermore, we are considering at most one scalar 
at one time and set the two-form B^ u = such that the non-abelian DBI action [43] 

S = -T p I STiV- det(P[E + E{Q~ 1 - t)E] + 2irl*F) det(QE) 

J Dp 



+ M p / TrlP 

I Dp V 



A2irliF 



reduces to 



S 



Dh 



STrV '-det(P[G] + 2ttP s F) 



(2.4) 
(2.5) 

(2.6) 



evaluated in the D3 background (2.1). 

Here, we want to consider the scalars in the z direction (i.e. along the D3 branes), 
<]? 2 , and the spherical direction $ sm£l (away from the D3s) that describe the position and 
separation of the probe branes. The Pauli matrices are (with an additional normalization 
factor) : 




n 




72 




73 




(2.7) 



with the properties 

[T a ,n] = ie c a bTc , {T a ,u} = <Wo [r a , t ] = and {r a ,r } = r a . (2.8) 

In particular, we want to consider configurations that are isotropic in the defect directions 
and static, such that the only coordinate dependence is on u and we also we want to turn 
on only one of the non-trivial SU{2) generators and as our results will be invariant under 
an overall SU(2) rotation, we can choose e.g. t\. Choosing the parametrization of the S 5 
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sphere as d£l\ = d9 2 + sin 2 OdQ?, + cos2 GdQ 2 anci putting the D5 on the second 2-sphere, 
we will choose the Ansatz 

<S> Z = 2^(u)ti and $ sine = 20(«)n . (2.9) 

This gives the induced metric (if we fallaciously evaluate the multiplication of the SU(2) 
generators already now) 



ds 2 = ^2To(-(l-u i )dP + d^+(^^+^ 2 + ^^)du 2 ) + L 2 (l-<p 2 )dn 2 , 



( l^ „A\A& i a~.1 i ( ^ i „/.'2 i U 2 (f)' 2 

^ J uu j -f- 1^ — qj )i 

(2.10) 

where (•)' denotes d u (-). Just as a motivation let us look at the case of the <£> z scalar only. 
The solution is just 

d u ip = C (2.11) 

V(l-U 4 )((l-U 4 )-C 2 u8) 

for some constant c and some U(2) valued matrix r. This tells us that without additional 
fields, the only embedding in which the branes fall into the horizon is the trivial one. In 
the only other case, the branes extend out to tp = ±oo at u A = + y / ' ^ + J. In order 
to find non-trivial embeddings, we should then turn on also a gauge field that does not 
commute with the scalar. 

To get a simple physical interpretation, we choose the gauge field to be proportional 
to t 3 with the Ansatz 

L 2 

A = — p (u) T3 dt , (2.12) 

7Tl s 

which can be interpreted as the gravity dual of the isospin density [37]. As the isospin 
number is essentially given by the appropriate number of string end/ "start" points on the 
probe branes, the isospin density and isospin chemical potential are then given by 

1 2 

hm d u p{u) = p = rjv Piso and limp(u) = p. = -j^zPiso ■ (2.13) 
mo 2iV c i z u~>Q \/\T 

Taking into account the newly arising commutators through the covariant derivative T>^& = 
d^<& + i[Au, <&], the pullback of the background metric (2.2) becomes 

ds 2 = 2^ (-{I - u 4 )dP + d5% + -^—] T + 
u z \ 1 — vr J 

L 2 

2^r (-2{iPpTz){il)pTz)dt 2 + Aii^pra^'njdtdu + 2(^ / T 1 )(iP'T 1 )du 2 ) + 
2L 2 

^—^ (-2^ P T 3 )(<PpT 3 )dP + 4i{<P P T 3 , (j>'Tl}dtdu + 2(cj)'T l ){<j)'T 1 )du 2 ) 

+ 2r L 2 (l - (j) 2 )dnl , (2.14) 

where we kept the generators explicitly because one has to take the symmetrized trace. 

To complete the picture, we will also alternatively to the scalars consider a "magnetic" 
field, taking the Ansatz 

L 2 

A = ^ {2p(u)r 3 dt + u(u)ndx) . (2.15) 
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In order to determine the phase diagram, we obtain the free energy in the usual way 
from the regularized euclidean on-shell action. By studying the variation of the fields, one 
can trivially see that this is actually a function of the chemical potential, rather than the 
density, so we get [44, 45, 46, 47, 48, 41]: 

n(2» = T/r 9 - . (2.16) 

The counter terms that we need to add are dictated to us by requiring consistency of the 
variational principle and give in analogy to the abelian case [49] 

hdy. = -\Vl+l<t> 2 V7 ■ (2-17) 

As we will see later, it turns out that the other scalar ip does not contribute to the boundary 
term. The Helmholtz free energy can then be obtained from a Legendre transformation, 

F(T,p) = Q + pn. (2.18) 

As it seems most physical, and also for reasons that we will see later, we want to fix the 
isospin density and hence work in the canonical ensemble using the Helmholtz free energy. 

To obtain some dimensionless quantity that naturally appears in terms of p or p, and 
the dimensionless boundary values for u), <j> and ip, we can define the usual dimensionless 

F f2 

F := — and := -= . (2.19) 

v^iVcT 3 ^/XN C T 3 

Since it is for most processes however more physical for a superconductor to fix the particle 
number rather than the temperature, the dimensionless quantity 



F 2J2N~ C 

(2.20) 



p 3 / 2 vV/2 

is the more physically relevant one. In oder to obtain a positive value, we define the energy 
gap between the unbroken and broken phases as 

— -fjmbrofcen ^broken • (2.21) 

In appendix A, we will also refer to the Gibbs free energy, where we define the energy gap 
exactly the other way round, Af2 = VL broken - ^unbroken- 

3. Perturbative expansion 

In the following, we look for convenience only at the leading a' expansion, where we scaled 
the scalars $ a such that they appear at the same order of a' as the gauge field A. The 
action for this pertubation is then 

5 = -^L^\N C J d 4 CTv^STr^2P M $ sine ^$ sine -4$ sine $ sin9 + ^P M $ z ^$ z + F 2 ^ 

(3.1) 
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where the metric g is just the induced metric without the scalar perturbations. Substituting 
in the Ansatz (2.9,2.15) and the metric (2.14), we obtain 

(3.2) 

The volume factor V2+1 for the dimensionless "flat" coordinates that we introduced here 
gives rise to a factor of 2-/r 3 T 3 in the densitized euclidean action. 

The solution in the case where we turn on only the isospin potential is trivially 

p(u) = q{l-u) , p = p = q (3.3) 

for some number q. This gives the free energy F = — \p 2 - 

Before we start using this action to obtain the phase diagram, let us review the appro- 
priate first basic examples in the literature. 

3.1 P-Wave superconductor 

The p-wave superconductor was found in a bottom- up approach in [17] and implemented 
in a D3-D7 system with 2 D7 probes in [34, 36, 35] and in our setup in [26]. Up to a 
difference of 1 — u 3 vs. 1 — u 4 in the blackhole factor the bottom-up and top-down models 
should be the same at the perturbative level. 

This case is called P-wave, because the symmetry breaking is caused by a vector field 
and the resulting spectral curves are highly anisotropic [17]. As an Ansatz for the gauge 
field one can then choose A = p(u)T%dt + uj(u)ridx, where oj(u) will play the role of the 
symmetry breaking condensate, and one could certainly choose any other combination of 
non-commuting generators. 

In [17], the equations of motion are (in our notation): 

(«V(«))' = ^P(«) and ((1 -«V(u)) = ^"(«) ■ (3-4) 
1 — u 6 1 — u 6 

The expansions near u = are 

p(u) = p — pu+ ... and u(u) = oOoo + u + . . . (3.5) 

and close to the horizon they are: 

p(u) = - q (u - 1) + . . . and uj(u) = ujq + 0(1 - u) 2 , (3.6) 

the first of which arises from requiting \A\ to be finite on the horizon, and the second 
because the equation of motion reduces to first order on the horizon. For a consistent field 
theory, uj^ must vanish. Hence we are left with a one-parameter solution that traces out 
a line in the (q, loo) plane or the (Q,£li) plane. Obviously p corresponds to a chemical 
potential and Q to a charge density, and oc< > makes a good candidate for the order 
parameter (see eq. (2.13)). In practice we will always represent the condensate scaled by 
an appropriate factor of the density, in this case C,uj = 



S 



VXN C V 2 
2vr 3 



+1 



du- 
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In our case, the equations of motion can be straightforwardly obtained from (3.1) and 
become 

(u 2 p'(u))' = y^pW) and ((i-«Vw) = r=^w ( 3 - 7 ) 

but the near- horizon and asymptotic expansions remain the same. 
3.2 S-Wave superconductor 

In the case of the s-wave superconductor, the symmetry-breaking operator is a charged 
scalar in a charged planar black hole background and it was first found in a bottom-up 
setup in [19, 20]. Later, for example various M-theory compactifications were studied as 
to whether they obey superconducting behavior in [33] . 
The action in [19] was 

C = F 2 + - |£>$| 2 (3.8) 

4 L 2 

with the equations of motion 

= 8 U ({u~ 2 - u)d u 0) + u~ 4 ( u _f_ - +2)<t>, (3-9) 

= d 2 u p - 2^-^ rp , (3.10) 

u*{u z — u) 

the asymptotic expansion 

p = p — pu + . . . and 4> = Ciii + (2U 2 (3-11) 
and the near-horizon solution 

p = - q ( u - 1) + . . . and = 0o + 0(1 - u) 2 . (3.12) 



Furthermore, consistency and stability apparently requires to set either = or (f>2 = 0. 
This gives rise to spontaneous symmetry breaking, if one sets e.g. Ci = an d considers (j 
at the same time as a source. 



3.2.1 Flat scalar 

As we discussed in section 2, we have two scalars at hand that we can use for the symmetry 
breaking, the one corresponding to the separation of the probe branes in the z direction 
along the D3 branes, and the one on the sphere separating the probe branes from the D3's, 
corresponding in the usual picture to a finite quark mass [8, 7, 47, 50, 51, 52, 45] or chiral 
symmetry breaking [53, 10, 9, 54]. For now, let us choose the former case in which the 
Lagrangian in the action (3.2) becomes 
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The equations of motion are then 

= d 2 u p- ui( f_ u4) P and (3.14) 

o 2 

= d u ((n" 4 - l)d u 4>) + • (3-15) 

1i 4 — 1 

Now there is a significant difference to the usual S-wave model, since there are two more 
powers of u~ l inside the second (gauge-covariant) derivative of cj). Consequently, the asymp- 
totic scalings are different and read now 

p = [i — pu + . . . and tp = c-p~ 5 / 2 ((3 + u 2 p 2 ) sin up — 3up cos up) + . . . = Cipu 5 + . . . 

(3.16) 

and we define = S$. In principle, there is one another solution, where p is approximately 
divergent, p ~ e^l u and if) is a decaying rapidly oscillating mode, 

ip ~ c+p~ 5 / 2 ((3 — u 2 p 2 ) cos up + 3up sin up) , (3-17) 

but we discard this solution as it should be non-normalizable. The near-horizon solutions 
are 

p = -q( u _i) + ... and ^ = ifj + 0(1 - u) 2 . (3.18) 



i. 



11=0, r/, -11.9 






n=l, 4 = 17.0 






n=2, 4=22.3 






n-3, 4 -27.7 






„=4, 9l =33.0 






n-5, 4-38.5 






n=6, 4=44.0 






n=7, 4=49.4 



















0. 0.2 0.4 0.6 0.8 u 1 



Figure 1: The solutions for il)(u) close to the 
critical temperature, normalized to 1. 



<*0 




0. 2.5 5. 7.5 10. 12.5 15. q 17.5 

Figure 2: The solution space (lines) of the 
internal parameters. The lines correspond 
to the different numbers of modes, as in 
the other figures. The values for q c are now 
{11.9, 17.0, 22.25, 27.65, 33.0, 38.5, 43.95, 49.4, 54.9, . . .} 



To find solutions, we first tune i/jq to a small value and study the linearized equations 
of motion, 

o = d u {(u~ A - 7i)c\v) + q ( ! ~ r ^ and P = tfC 1 - «) ( 3 - 19 ) 

' u 4 — 1 

and tune q such that we satisfy the asymptotic boundary condition. We find that there are 
slowly oscillating solutions for tp with n £ Z^J" modes at various critical q c that we show in 
fig. 1. 
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In terms of the internal parameters q and V'O; we find an interesting attractor behav- 
ior that we show for in fig. 2, as they all flow in an in-spiraling behavior to (q, ?/>o) ~ 
(6.409100,0.8294795). It seems that every valid solution, at least in the first quadrant of 
this phase space, that matches the boundary condition (3.16) will flow to this point in the 
limit of small temperatures, and because of the symmetry of the action, every quadrant 
will have such an attractor point. 




Figure 3: The value of the condensate for the Figure 4: The gap in the free energy between 
odd-mode solutions are negative, so the abso- the unbroken and superconducting phases, 
lute value is shown here. The connected line 
corresponds to values obtained at u ~ 0.1/i~ 4//5 
and the dots to«~ 0.05^ -4 /' 5 . 

Obtaining the value of the condensate C/, is somewhat non-trivial as this quintic mode 
is hard to detect. At very small values of u, the solution is dominated by a small numer- 
ically remnant constant term, and the range over which the quintic term is dominant is 
very small as the solution quickly levels off into the "finite-ti" regime. The most reliable 
method is to fit both asymptotic modes up to u ~ 0.1/2~ 4//5 . Still, there is some numeric 
noise, and the results in the and 1-mode solutions have already ~ 2% systematic error. 
To suppress the noise, we consider only values for which the coefficient of the constant term 
is sufficiently small. The results illustrating this are shown in fig. 3. The result for the 
1-mode solution remains below the 0-mode solution independent of systematic errors. An 
interesting observation is that we had to plot T 5 / 2 ^, to show a finite value at small tem- 
peratures, so the condensate diverges proportional to T~ 5 / 2 . This however is no surprise, 
since we are considering an incomplete action that can give rise to singular solutions, and 
this behavior should disappear when one considers the full DBI action. 

Finally, we look at the energy gap in the free energy in fig. 4. We find that we have 
indeed a second order phase transition with a continuous derivative of the free energy at 
the onset of the symmetry breaking solutions, and the 0-mode solution is the preferred one. 
There are two puzzling aspects though. Firstly, we have plotted TAF, rather than just F 
which means that at small temperatures, the gap in the free energy is actually diverging, 
which is clearly not physical. Also this divergence should disappear if one considers the 
full DBI action. Secondly, the value of the energy gap is approaching the same value as 
T — > 0. Maybe this is related to the attractor behavior, or it may also be some restoration 
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of supersymmetry at vanishing temperature. 



3.2.2 Compact Scalar 

Now, let us consider the second scalar $ sin9 ; corresponding to the separation of the probe 
branes from the D3 branes. Taking hence only the fields <f>(u) and p(u) in the action 3.2 
to be non-vanishing, the relevant lagrangian is 

Hence, the equations of motion are: 

= d 2 uP - u2( f_ u4) P and (3.21) 



2 



= d u ((u- 2 - u 2 )d u <f) + P J > + . (3.22) 
' u z (u 4 — 1) 

In their structure they are very similar to the equations (3.9) of the usual holographic su- 
perconductors [19, 20] and differ just by the factors and by the power of u in the blackening 
factor - the latter since we are considering the induced metric on the probe brane in an 
AdS§ black hole background rather than a planar AdS^ black hole. 
The asymptotic scaling is now just as in [19, 20] 

p = p — pu + . . . and 4> = (\u + C2U 2 (3.23) 

and the near-horizon expansions are 

p = -q(u - 1) + ... and = O U + \{u - 1)) + 0(1 - u) 2 . (3.24) 

Now there are two possibilities to consider either the operator corresponding to (j, as a 
source, or the one corresponding to £2 an d require them to vanish accordingly. Again, we 
can find the critical values of q by solving the equations of motion for small values of (f>, 

= d u ((n~ 2 - u 2 )d u 4>) + g2 2 ( / 4 ~ M f> + , p = g(l - u) , (3.25) 

u z (n 4 — 1) 

for modes with d u (j)(u)\ u ^ = or d 2 (f)(u)\ u ^ = 0, respectively. The scalar field is then 
just a decaying oscillating mode as shown in fig. 5. 

Then, we can increase the value of the scalar on the horizon, (fro, and track the corre- 
sponding value of q c , as shown in fig. 6. Now, find again a kind of attractor behavior with 
a divergence (fio ~ —0.33 In q + const., where the constant increases with increasing modes. 

Looking at the result in terms of the physical parameters in figs. 9 and 7, the con- 
densate Ci = ^ or ^2 = p and the temperature T, we find a picture that is qualitatively 
similar to the previous case - but with much larger critical temperatures. Now, only the 
th mode of the (2 = case diverges, proportional to T 1 / 3 . Interestingly, (2 flows precisely 
to ±1 for all modes at small temperatures. 
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0. 5. 10. 15. 20. 25. q 

Figure 6: The internal parameters </>o and q c for the first few mode solutions. The values for Ci = 
are shown dashed. The inset shows a logarithmic-linear plot. 




Figure 7: The value of the condensate £2 as a Figure 8: The superconducting energy gap for 
function of temperature. the first few modes. 



Finally, we can again look at the energy gap in fig. 8, where we see clearly again 
a second order phase transition, and the fact that the £2 = 0, n = mode is the ther- 
modynamically preferred phase. We find again the divergence of the energy gap at small 
temperatures, and the flow to F = 1/(2T) at small temperatures. 
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Figure 9: The value of the condensate £1 as a function of temperature as a linear and double 
logarithmic plot. Again, the value of the condensate is negative for odd n and only the absolute 
value is shown. 



Now, let us compare the different super- 
conducting phases to determine which one 
is the stable one. Using again the canon- 
ical ensemble, we see in fig. 10, that the 
dominant condensate is always the s-wave 
£i- For comparison, we also show the n = 
mode of the p-wave, and it turns out that 
this is only metastable and should decay to 
the s-wave. 

Interestingly, the p-wave and the £2 s- 
wave condensates behave very similarly. 

Certainly, we assumed here that the canon- 
ical ensemble is relevant as fixing the parti- 
cle number is a sensible thing to do for superconductors. For completeness, we also provide 
some results related to the grand canonical ensemble that may be relevant in QCD contexts 
in app. A.l. 

4. Full DBI Action 

Now, let us consider the "full" DBI Action. Certainly the non-abelian DBI Action (2.4) is 
not complete and receives corrections at higher order in a' [55, 56, 57, 58]. In any event, 
so far it is not possible to compute even the non-abelian DBI action exactly because of the 
symmetrized trace. The symmetrized trace prescription requires us to take a symmetrized 
average over all orderings of T>^<&, F^ v (and generically [<fr a , $ ]). Hence, in oder to obtain 
an exact solution, we would need to do some resummation, which will probably be possible 
only in the simplest cases. In order to still preserve the features of the regularity of the 
DBI action and the geometry of the S 5 sphere, we will do the following: 



3.3 Comparing the Condensates 




Figure 10: Energy gap comparing the different 
condensates. 
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Defining the symmetrization operator S, we can write the DBI action schematically as 



S = - STtV- = - / TrSf- = - / TtSJS{-), (4.1 



where in the second identity we used S(S(f)S(g)) = S(fg) and assumed that we can 
expand the square root. Even inside the square root, one can go as far as symmetrizing 
any pair of the "atomic" expressions. The "approximation" that we will do is to evaluate 
this "first" symmetrization in the matrix representation (2.7), i.e. we will replace any pair 
of such objects with their anti-commutator. Since the action (2.6) contains these terms 
only at even order, the expression inside the square root will even be proportional to the 
identity, and one can readily compute the result. 

That this procedure is not exact can be seen at the following example: S{t\T2T\T2) = 
S(t\TiT2T2) because of the symmetrization, and we still have S{S{t\T2)S{t\T2)) = S (S (tiTi) S (T2T2)) ■ 
If we fallaciously evaluate the anticommutators in each side using the representation (2.7), 
the left hand side will vanish, but the right hand side will be proportional to the identity. 
As our "approximation" is however identical to the so-called symmetrized trace prescrip- 
tion that has been used in the literature [35] and it preserves the features of regularity and 
the spherical geometry, we will use it for the rest of the paper. Because of the commutator 
term [<3? a , <3? 6 ] in (2.4), let us not write down the action now, but postpone this to the 
separate cases. The induced metric simplifies in this "approximation" to 

ds 2 = 2^t (-({I- u 4 ) + P V + u 2 /^] dP + dP 2 



U A \ \ 1 — Z 

+ + ^ + ^2) ^ 2 ) + 2r L 2 (l - ^)dtf 2 . (4.2) 

Before we start finding the solutions in the broken phases, let us find the solution of 
the unbroken phase. In this case, the action in the parametrization of (2.12) becomes after 
an integration over the sphere 

s = -^yftNcJ d^V^V • (4-3) 

Note that this is a somewhat unusual normalization of the gauge field, arising from the 
definition in (2.12). The solution for p is 

p' = — / = = = j P = %\/iq sinh -1 \/iq, i) — F(i sinh" 1 \fiqu, i)\ , (4.4) 



where the parameter q is related to the isospin density through p = 2q as there are now 
no sources. F(-, ■) is the elliptic integral of the first kind. At large values of p, this gives 
us P ~ .'f- I 1 (z) 2 ( ? 3//2 ; whereas at small values, we have p ~ p, as expected. The value of 



the action becomes 



S = -^VXN c V 2+ i x ? ^^^(isinh- 1 y/iq,i) - ^1 - g 2 ) , (4.5) 
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where we included the usual counter term 5^/7 and introduced a volume factor V2+1 for the 
dimensionless "flat" coordinates, that gives rise to a factor of 2-7r 3 T 3 in the densitized eu- 
clidean action. In the limit of large q, this expression evaluates to S ~ — ^^^ 2+1 (|) 2 

and at small q, it becomes just the usual conformal S = — v ^^f 2+1 §■ 
4.1 P-Wave 

Now, let us consider only the "magnetic" field uj{u)T\dx as in section 3.1. The action 
becomes now in the parametrization of (2.15) and after an integration over the sphere 

v^c r a~„ 1 A, ; .4 ^ ^ o ^ .. 4 p 2 ^ 2 



5 = d*a- A] j4 + uH(l-u^-^)-u*f-^, (4.6) 

such that equations of motion for p(u) and uj(u) become 

d u , P ' ^ (4 

J4 + u 4 ((1 - u 4 V 2 - p' 2 ) - u 4 ^r (l-n 4 )^/4 + u 4 ((l-u 4 )u; 2 - p' 2 ) - 

M , "' (1 -" 4) J ~ , ^ ..-(4 



4 + u 4 ((l-u 4 V -p /2 ) -n 4 ^/ (l-n 4 )^/4 + n 4 ((l-u 4 )a; 2 - p /2 ) - 



2,.,2 
4 



Near the horizon, the boundary conditions can be found by an expansion of the equations 
of motion around u = 1 as 

P = - ? P=(l-u) + ... and u = uo- U ^{ n =^\{l-uf + ... , (4.9) 

where as usual the vanishing of the "electric potential" on the horizon is required by 
regularity. We chose the parametrization of the electric field analogous to the solution 
(4.4) for the gauge field in the unbroken phase in order to ensure regular solutions and also 
for practical purposes for the fine-tuning near , q = 1. Also physically, q still represents 

the flux through the horizon, or the isospin density that is not carried by the condensate. 
The asymptotic expansions are the same as in the perturbative case (3.5), p = p — pu + . . . 
and lo = ( LU u + . . .. 

In fig. 11, we look at the relation of the internal parameters ujq and q for the first three 
modes of the spontaneous symmetry breaking solutions. As usual, the value of q decreases 
with an increasing horizon value of the symmetry breaking field. Now however, because of 
the "shape" of the DBI action, the critical values of q are separated by a large factor rather 
than just equally-spaced. The critical temperatures are correspondingly suppressed, at 
values f c ~ {0.39327, 0.094872, 0.019068, . . .}. The divergence of u at small temperatures 
is logarithmic with ojq ~ — 0.87 In q + const, or expressed in another way, the value of q is 
exponentially suppressed at large ojq. This just reflects the fact that most of the density is 
carried by the condensate in that limit. 

The profiles of uj(u) and p(u) are shown in fig. 12 for some selected values of the 0-mode 
solution. We find the usual "positive feedback" between the symmetry breaking field and 
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Figure 11: The relation of the internal parameters for the first few modes of the spontaneously 
symmetry-broken solutions for the p-wave condensate. 




Figure 12: The n — mode solution of the p-wave superconductor for different choices for the 
boundary condition: Symmetry breaking field (left) and isospin potential (right). The black narrow 
line shows the isospin potential for q = 3.14 for reference. Curves for larger values of q will be below 
this curve. 



the isospin potential. This also illustrates how the isospin chemical potential or density is 
supported by the symmetry breaking field: While the profile for the isospin potential for 
large values of q and small values for loq is similar to the one in the unbroken phase (as it 
should), for small values of q and large values of loq the potential (and its "electric field") 
are almost entirely sourced by the charged field, which in turn is approximately constant 
up to large radii (or small values of u). 

In fig. 13, we show the (absolute) value of the symmetry breaking condensate £ w = -7= 
in the broken phase. While the critical temperatures for the n th mode scale approximately 
proportional to 4 _n , the small-temperature limit of the condensate rather saturates - both 
are due to the form of the DBI action. 

In contrast to this behavior, the curves of the energy gap in fig. 14 are approximately 
invariant under an approximate scaling with the critical temperature. Also in contrast to 
the a' 2 expansion, the energy gap is finite at small temperatures. This is a clear feature of 
the regularity of the DBI action. 
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Figure 13: Absolute value of the supercon- 
ducting condensate for the first few modes of 
the p-wave superconductor. 



Figure 14: Gap in the free energy with re- 
spect to the unbroken phase. 



4.2 Compact Scalar 

Next, let us consider the s-wave condensate from the "compact" scalar that describes the 
separation of the D5 probes on the 5 5 . In contrast to the case of the a' expansion, this 
will contain some more interesting physics, as we have now black hole and non-blackhole 
embeddings in which the probe branes do not extend down to the horizon, just as in the 
usual massive probe brane embeddings (see e.g. [47, 48, 46, 41]). This may give us an 
additional first oder phase transition. 

In the parametrization of (2.9), the action becomes now 



2vr 3 J u A ^ ; \j V (l-<p 2 )(l-u 4 ) J V y l-<£ 2 
which gives us the equations of motion 

P '(l - 4> 2 ) 



(4.10) 



U 



4 1 - " 2 ( i-/)ff-u-) ) 1 + « a (i " " P' 2u " 



4p0 2 ((l-^ 2 ) + u 2 (l-uV 2 ) 



U 



U 2 {1 _ U 4 )(1 _ 02)^4 _ U 2_J^_^ ^ ~ U 2 {1 _ U 4)_^) _ p ,2 u 4 

4((l-^)(l-« 4 )-«W)^' 



(4.11) 



4 u V^ + 2(l-</) 2 ) +4u 2 <P' 2 (l-n 4 + «V)-«V 2 (l-^ 2 ) 



■ (4-12) 



«*(! - ^) A /4 (l - « 2 (i3#)tW) ) I 1 + « 2 U " « 4 ) A- ) " ^ 



n'2,,4 



The asymptotic embeddings are again p = fl — pu + . . . and </> = Ciu + C2W 2 + • • •, where 
either £1 or £2 is considered as a source and hence required to vanish, and the other one 
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is considered to be the condensate. Again, we define (i = and £2 = p- For the "IR" 
boundary condition, we have to consider the black hole embeddings and the non-blackhole 
embeddings separately. 

4.2.1 Black Hole Embeddings 

The usual near-horizon expansion around u = 1 gives us 

' = 7TT? (1 -«) + - * = *„ (l + + . . . , (4.13) 

where we chose again the parametrization in terms of q to obtain a regular solution and 
reflect the flux sourced by the black hole. 
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Figure 15: The horizon value of the scalar field as a function of the "isospin flux' through the 
horizon as a linear and a linear- logarithmic plot for the n — and n = 1 modes. 

Looking at the "IR" parameters for the n = and n = 1 mode solutions in fig. 
15, we find again that the ordering of the critical values for q is the same, and they are 
exponentially separated, rather than linearly as at order a 12 . There is an attractor at 
q = and 4>q = 0, which just reflects the geometry of the sphere. It tells us that all of 
the charge density becomes supported by the condensate as the embedding approaches the 
critical embedding in which the brane touches the horizon at only one point. Intuitively, 
this arises because an embedding with finite flux out of the horizon would give a singular 
value for the action at this point and is hence avoided by the DBI action. 

4.2.2 Minkowski Embeddings 

Now, let us look at the case where the branes do not extend to the horizon, but the black 
hole is located "in between" the branes. In order to find the boundary condition at the 
point uq at which the size of the S 2 sphere vanishes, — ?■ 1, we substitute the Ansatz 

<j> = 1 - (p' (uo-u a ) and p = q(u - u) 13 (4.14) 

into the equations of motion (4.11). In any case we should assume o, b > in order to obtain 
a sensible solution with the given Ansatz, and furthermore we will assume the limitation 
of all the possible permutations of 



{2/3 - 2, a - 1, 0} < {2 - a, 4/3 - 2 - a, 2a - 1, 2/3 - 1} 



(4-15) 
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Note that the solution for the abelian case, 6 = 1/2 and a = 1 violates that bound 
marginally. Under the condition (4.15) we get now the four leading orders 

= u%q 2 p(u - n) 2/3 " 2 (2a 2 - a/3 - 2f3 2 ) - 4a 2 ^ (n - u) 2f} ~ 2 (3 - uf) + (8 + 4a)(4aifi) 
= (n - u)' 3+a ~ 2 (2a 2 - a/3 - 2/3 2 ) . (4.17) 



The non-trivial solutions to these equations are /3 = 1, (j)' = 3u ^_ u 5 and a = — 1(1 db y/V7), 
such that the only solution that also satisfies a > (and the bound (4.15)) is 

<fi ~ 1 k( u o ~ u ) an d P ~ q (i*o — u ) 2 • (4-18) 

3n - ng 

It was argued in the abelian case of the D3-D7 system in [52] (and for the D3-D5 system 
in [41]) that the D3-Dp strings that are necessary to source the flux in the Minkowski 
embedding of the probe branes would be in disbalance with the tension of the Dp-branes 
at no and cause the Minkowski solution to collapse into a black hole solution. However, it 
turns out that the electric flux 



/ a / / A \ / -t i 0\ O (ON /-. TO i lO / -, A \ \ tO /-, A\ V / 



68 u A t y/4 ((1 _ u 4)(! _ ^2) _ p 202) (1 _ 02 + ^2(j _ „4)) _ ^ _ u 4 

does vanish in our case near n ~ no as 



6C V17 



<55 U A 4 4n 0v /l 




17+1 



oc q- -== f W 5 .x(no-n) 4 . (4.20) 



This just means that the flux on the probe brane is entirely sourced by the scalar conden- 
sate. Hence, the parameter q is now related to the density of Dp-Dp strings near no rather 
than the number D3-Dp strings in the abelian case. In fact, one can consider the variation 
of the flux, d u sd^At t° describe the density of the source-Dp-Dp strings, and from (4.20) 



V -L ' — o /l V Q 

we see that it scales as q(no — n) 4, . As v 4 ~ ~ 0.28 and hence the density vanishes 
as n — > no, this is consistent with the fact that the scaling of the probe brane embedding, 
(4.18), does not depend on the value of q. 

To discuss the dependence of q on no, we have to keep in mind that for large no close 
to 1, the brane profile has some non-trivial "trumpet- like" shape (i.e. it extends almost 
down to the black hole), and for small values of no, it spans almost "straight" from no to 
(i.e. it is approximately flat, at large radii). This tells us precisely why the value of q that 
we need to obtain the 0-mode solution for </>2 = vanishes as no — > as shown in fig. 16: 
The attraction from the black hole is small, such that the boundary condition is almost 
matched by the geometric part of the DBI action and one needs little flux to enforce the 
boundary condition. Some more consequences of this will appear later. However for the 
higher mode solutions and for (f>i = 0, q becomes large as we need a large flux to support 
the profile, also because the "space" that is available in the inverse radial coordinate n to 
"bend' the brane vanishes. It is somewhat surprising though that the scaling exponents 
around no — > seem to be non-rational. Because of our definition of q that does not contain 
the other factors in (4.20), q vanishes oc y/l — Uq as no — ^ 1. 
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Figure 16: The horizon value of the scalar field as a function of the Dp-Dp string density parameter 
q as a linear and a double-logarithmic plot for the n — and n, = 1 modes in the Minkowski 
embedding. 



While in the black hole embedding phase the isospin density is also carried an appro- 
priate number of free quarks and anti-quarks that are in equilibrium with "mesons" , now 
all the isospin density is carried by mesons. 

4.2.3 Combined Results 




Figure 17: The 0-mode profiles of the probe brane embeddings for both Minkowski and black hole 
solutions for ( 2 = (left) and Ci = (right). 



In fig. 17, we show the brane profiles of the 0-mode solutions for the two different 
condensates for a few values of 0o and uq , respectively. We see that the profiles for 0o and 
uq close to 1 actually cross close to the horizon, and also the curves for (J)q ~ 0.5 and (fio ~ 1 
are similar in the UV (small u) . These features are reminiscent of the meson melting phase 
transition [47, 41] and its many aspects/relatives in different regimes (see e.g. [10, 9]). 

Just as in the case of the a' 2 expansion in section 3, the n = mode is again different 
from the other modes, as we show in fig. 18. Now, it starts off again proportional to 
y/T c — T, and there follows a multivalued behavior that is typical for black hole - Minkowski 
transitions. However rather than extending to T = 0, the phase "ends" at some finite 
temperature T* at which the value of the condensate diverges, approximately proportional 
to (T* — T) -1 / 4 as we show on the right plot. Prom the gravity side, this is not surprising, 
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Figure 18: The value of the condensate £i for the first two modes. The insets on the left show 
details of the n = 1 mode and the crossover between the black hole and Minkowski embedding. 
The plot on the right shows the condensate scaled by (T* - T) 1 / 4 for some T* = 0.905956 that is 
indicated by the black line. 



as at vanishing no, Ci naturally diverges as Ci ~ l/«o- Since the scale that sets the 
temperature is just given by the isospin density, i.e. by the "electric flux", and we only 
need some finite flux to support this embedding as discussed above, the phase ends at this 
finite T*. Usually the phases extend to T = 0, where the divergence of the condensate 
is canceled by its temperature scaling due the to dimensional factor. We will discuss the 
relevance of this effect further below. 

The n = 1 mode starts rather unusually proportional to y/T — T c - probably a sat- 
uration effect due to the regularity of the DBI action. Then it has the transition to the 
mesonic phase, which extends to T = 0. The fact that the value of the mesonic phase 
is almost constant is qualitatively not surprising, as in the mesonic phase most the free 
quarks have already condensated. While it is certainly still not impossible that there is 
a smooth density of D3-Dp strings that can potentially condensate, these stretch very far 
and are hence energetically disfavored 
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Figure 19: The value of the condensate £ 2 as a 
function of temperature. For the n = 1 phase, 
we could only find the black hole embedding. 



Figure 20: The gap in the free energy for 
the first few modes of the "compact" s-wave 
solution. The inset shows the energy gap scaled 
by (T* - T) 3 / 4 . 
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For if2 = 0, we could find for the n = 1 mode only the black hole embedding, as the 
Minkowski embedding was beyond reach of the numerics for most values of uq. The results 
are shown in fig. 19, and we find that now already the n = mode is similar to the n = 1 
mode of the other scalar. This is because one needs a large charge density to "bend" the 
brane profile, and this is also sourced in regimes away from the UV. Overall, we also see 
that the results get somewhat noisy, especially close to the critical embedding at no = 1 or 

00 = 1- 

Looking at the energy gap for the Helmholtz free energy in the canonical ensemble in 
fig. 20, we see again that the (f2 = 0, n = phase does not extend out to T = 0, but it 
diverges at T* . The power is not quite clear and can range from —0.66 to —0.75, but if 
one wants to have the same T* as in the diverging condensate, values close to —0.75 seem 
to give better fits. The "drop" that we see in the right plot in fig. 20 can be due to a 
subleading term that is proportional to (T* — T) 1 / 4 . The ratio of the divergences of the 
energy gap and the condensate is then 3, consistent with the ratio of their dimensional 
scaling. The geometrical reason is the large volume factor of the probe branes that extend 
only to some small uq (i.e. large minimum AdS radius), so it is perfectly consistent, and 
the problem again only arises because the flux needed to support this embedding is finite, 
rather than diverging. This divergence seems to be an artifact of the adapted symmetrized 
trace prescription, as we will discuss in the conclusions in section 5. Before the divergence, 
we see a second order phase transition and some multivaluedness, which is in a similar form 
also seen in the usual black hole - Minkowski transitions. Rigorously speaking, the first 
turning point would give rise to a th oder phase transition, which would not appear if the 
Minkowski phase would extend to vanishing temperatures as this part of the black hole 
would not be realized and we would have a first order phase transition to the Minkowski 
phase before. 

In the higher modes, the superconducting phase transitions "reverse", such that the 
black hole embeddings are always metastable and we rather have a direct first order phase 
transition to the Minkowski phase. The change of the order of the phase transition is 
qualitatively not surprising, as in the black hole phase there is an equilibrium of mesons 
and free quarks, so the condensate can condensate smoothly, but in the mesonic phase, 
most of the free quarks have already condensated, so one cannot have a continuous phase 
transition. 

4.3 Flat Scalar 

Finally, we consider the case of the scalar & z that originates from the embedding of the 
probe branes in the flat (field theory) directions. Using the parametrization (2.9), the 
action is 
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which gives us the equations of motion 

(4 - (l + V^' 2 ^) - «V 2 4(1 - (4 - g£) (l + #21-1) _ u y 2 

g u ^(4(l-n 4 )-pV) = -/^ 2 (4 + y/ 2 (l-n 4 )) 

' (l+f^f)-«V 2 4(1 - U 4) J(4 - g£) (l + # 2 ^) - «V* 

Again, we expand around u = 1 to obtain the near-horizon limit 

P = ^p=(l-u) + ... and ^ = Vo- ^(-fM=)(l-u)* + ... (4.24) 

and the asymptotic expansions are identical to the ones in the perturbative expansion, 
p = p, — pu + . . . and ip = Ctp u5 + • • •> that we discussed in detail in section 3.2.1. Also, we 
define again C-0 = %■ As in section 3.2.1, there are some technical problems due to this 
highly suppressed scaling of tp in the UV, so in particular the results for the n = 1 mode 
will be very noisy and only available in some regimes at all. This noise is why we will show 
the individual data points in some of the plots. 
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Figure 21: The boundary condition for the scalar in terms of the flux through the horizon. The 
right plot shows the details of the attractor region at small q. 

First, let us look at the "internal" parameters, the flux out of the horizon q and the 
"separation" of the branes on the horizon, t^Q in fig. 21. We find again an attractor 
behavior with a fixed point around (q,4>o) ~ (11.2,0.385), but now the "spiral" is heavily 
distorted and the critical value of q for the n = mode moved to a very large value of 
q c ~ 1423. For the n = 1 mode the numerics could not reach far enough to obtain q c . 

In fig. 22, we show on the left the profiles of tp of the 0-mode solution for choices of 
q and V'o near the minima and maxima of the attractor spiral. We see that the curves 
seem to approach some non-singular limiting profile and the large value of the condensate 
is obtained by following that profile more closely and shifting the oc it 1 behavior further 
into the UV. The temperature dependence of the condensate is surprisingly similar to the 
higher modes of the other scalar in figures 18 and 19, even though they are geometrically 
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Figure 22: Left: Embedding profiles for the "flat" scalar for the n = mode for several values of 
tjjfj and q near the maxima and minima of the attractor spiral. Right: Scalar condensate C/> as a 
function of temperature. For the n — 1 mode only a small regime is shown due to numerical issues. 

very different and also the internal parameters behave differently. This may be a universal 
feature of scalars in the DBI action. Physically this may arise because we are looking 
at very small temperatures, i.e. large dimensionless densities p, while the dimensionless 
flux parameter at the horizon, q, that describes the un-condensated isospin density quickly 
moves towards its finite attractor value. Hence as soon as the internal parameters are in the 
attractor region, most of the density has already condensated and chance the condensate 
is approximately constant. For the n = 1 mode, we show only a small region in this plot, 
as the numerically obtained value of (p§ becomes extremely noisy. 

As with the value of the condensate, 
also the energy gap, shown in fig. 23, is 
very similar to the other scalar, albeit with 
a much lower critical temperature and corre- 
spondingly lower value. Again, we see that 
the second order phase transition "inverts" 
and we have now a first order phase tran- 
sition directly to the regime in which the 
density has almost completely condensated. 
The curve for the n = 1 mode is again very 

Figure 23: Energy gap for the scalar condensate similar to the n = curve ' u P on a scalin S 
corresponding to the position of the probe branes with the critical temperature, 
in the field theory directions. 

4.4 Comparing the Phases 

In fig. 24, we compare the energy gap of the different condensates. We see first of all that 
the ordering of the critical temperatures does not change compared to the a' 2 expansion, 
but especially the smaller critical temperatures are further suppressed. In addition to the 
changes discussed already in the previous subsections, we see now that there seems to 
be a first order phase transition between the p-wave and £2 s-wave phases. Obviously we 
studied the different phases only separately, and if one considers the full system there would 
be some mixing between the phases. Also it is not quite sure whether such detail would 
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Figure 24: The energy gap 
of the different condensates 
using the "full" DBI action. 



persist if one considers the rigorous symmetrized trace and the higher order corrections of 
[55, 56, 57, 58]. This would probably also cause the £i s-wave phase to extend to T = 0. 

Further, we also see that while the T = value of the energy gap is now different for 
the different phases, the first derivative, and hence the amount by which the entropy is 
lowered by the symmetry breaking becomes 



AS 



dAF 
df 



(4.25) 



at small temperatures. This is not contradicting the third law of thermodynamics or the 
positivity of entropy, as it is well know that these holographic systems describe at small 
temperatures a quantum liquid with non- vanishing ground state degeneracy [59, 42, 60, 61]. 
In fact, we can write down the free energy in the unbroken phase, 



unbroken 



T I Jl-(2T 2 y 2 + 



<J>2 



i/(2T 2 )F(isinh- i Ji/(2T 2 ),i 



(4.26) 



which gives us S — > 1 as T — > 0. Hence the symmetry breaking precisely cancels the zero- 
temperature entropy, such that the ground state degeneracy of the broken phase vanishes. 
Inspecting the curves further for the heat capacity c v = df(F + TS) = Td^F, we see that 
in the cases of first order phase transitions the heat capacity is increased whereas for the 
second oder phase transitions it is lowered. 

As in the perturbative case, we show some results related to the grand canonical 
ensemble that may be relevant in some other context in appendix A. 2. 



5. Discussion and Conclusions 



In this paper, we studied the phase diagram of a system of two probe "defect" D5 branes 
in the presence of a non-trivial SU(2) world-volume gauge potential in an AdS§ x 5 5 back- 
ground that is created by a stack of N c 3> 1 D3 branes. The field theory dual is a 3 + 1 
dimensional SU(N C ), N = 4 SYM theory coupled to two flavors of charged fermions that 
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are confined to a 2 + 1 dimensional "defect" . Certainly, the system is in any case still invari- 
ant under an overall SU (2) f rotation, but one way to interpret the system is by choosing 
p{u)T^dt (with an inverse radial coordinate u) as a potential which corresponds to a finite 
isospin density. Having thus broken the SU(2) explicitly to a U(l), we considered three 
different "channels" of spontaneous symmetry breaking of this U(l) using fields that do not 
commute with that potential (e.g. proportional to t%). The candidates were a magnetic 
component of the world- volume gauge field, ui{u)T\dx, giving rise to a p-wave superconduc- 
tor, and two different scalars giving rise to s-wave superconductors: A scalar corresponding 
to the embedding of the probe branes in the flat (3+1) field theory directions orthogonal 
to the (2+1) defect directions but along the D3s, and one corresponding to the embedding 
in the "internal" S 5 , away from the D3 branes. In some sense they could be imagined as a 
relative separation of the probe branes (if we had chosen the scalar proportional to T3), but 
it should be rather considered as some entanglement at two relatively separated positions. 

Technically, we used the non-abelian DBI action that governs the world-volume fields 
on the probe branes in the probe approximation Nf <C N c and studied the symmetry 
breaking fields separately. As it is physically most sensible for our considerations of super- 
conductors, we chose to work in the the canonical ensemble, where we fixed the density. 
For completeness, we referred to what would happen in the grand canonical ensemble, that 
would be relevant e.g. in a semiconductor in app. A. 

Qualitatively, our setup could be motivated from condensed matter physics by con- 
sidering a 2+1 dimensional system with (at least) two separate species of particles that 
interact in the effective field theory - e.g. a multi-layered system like cuprate supercon- 
ductors or multi-layered graphene. Depending on which kind of potential we choose, the 
density that we turn on corresponds then to some non-trivial combination of the particle 
number of either species, not only the "isospin" density of particles in one species and 
antiparticles in the other. In the string theory side, the symmetry breaking condensates 
correspond to some combination of D5-D5 strings, which translate on the field theory side 
into bound states of particles of different species. The unbroken density however corre- 
sponds to D3-D5 strings, which just translates into appropriate numbers of free particles 
in the field theory. 

There may be some criticism over the name "superconductivity" , as on the field theory 
side the U(l) that is broken is only a global symmetry. First of all, the local structure 
is not relevant for the symmetry breaking pattern, but it is only important once one 
considers further perturbations or external fields. Also, there are many AdS/CFT papers 
that observed e.g. the delta function in the conductivity or the Meissner effect - which 
are characteristic for local symmetry breaking. This should be also straightforward to 
implement here, as we could add a term like A magn , = Br^xdy that will act like an 
external magnetic field and one can straightforwardly study the conductivity using linear 
response theory. For this, we would expect to obtain results similar to the ones in the 
literature. 

5.1 Perturbative Results 

First, in section 3 we considered an expansion of the non-abelian DBI action to the level 
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a' 2 in order to give us a quadratic action (with quartic coupling) of the relevant fields. At 
this level, there is actually no coupling in between the different condensates such that our 
approach of considering them separately is exact, and one could in principle construct a 
more complex phase diagram involving combinations of the different condensates though 
linear combination. Comparing all the energy gaps in section 3.3, we found that the 
different phases look very similar, up to a scaling of the temperature dependence, and 
they all flow to the same T = diverging behavior with no crossing of phases. This 
might be due on the one hand to some supersymmetry effect, on the other hand, as we 
saw a AF oc \/T divergence, this is more likely just an artifact of the incompleteness and 
irregularity of the a' 2 expansion at large fields. All these symmetry breaking phases clearly 
showed a "superconducting" second order phase transition. We also looked at the grand 
canonical ensemble in app. A.l, where we found that the phases with broken symmetry 
would be actually disfavored if we had rather fixed the chemical potential. All of the 
phases had various higher mode solutions, and for the "compact" s-wave, there were also 
two different boundary conditions that we could impose, depending on what we consider as 
the (vanishing) source and what as the (spontaneously appearing) condensate. The higher 
modes always had smaller critical temperature and a smaller energy gap. 

The condensate with the largest critical temperature and largest energy gap was the 
compact scalar in which the condensate corresponds to some asymptotic ri-valued "sep- 
aration" of the probe branes from each other and from the D3 branes, discussed in sec. 
3.2.2. This is similar to an effective mass that is generated for the fundamental fields but 
not necessarily for their bound states. This particular condensate (only the 0-mode) had 
a diverging value of the condensate as T — > 0, so the T = values of the condensates of 
the different modes were approximately proportional to ra -1 / 2 . The condensate dual to the 
other boundary condition had a much lower critical temperature and showed an interesting 
flow of the value of the condensate to the same T = value for all the different modes. 

The p-wave condensate of sec. 3.1 had the second highest critical temperature, but 
we did not discuss it in detail as it has been readily studied in the literature [26]. The 
s-wave phase corresponding to the separation in the "flat" directions that we discussed in 
section 3.2.1 had the lowest critical temperature. While the "internal parameters", the 
horizon value of the "isospin electric flux" , i.e. the uncondensated part of the density, and 
the horizon value of the symmetry breaking field had a logarithmic dependence at small 
temperatures in the other cases, the scalar of the flat direction was special. There, we 
found an interesting attractor behavior to a finite fixed point. The condensate value of this 
scalar remains finite as T — > - computing it was numerically non-trivial due to the high 
order of the asymptotic scaling. 

5.2 Non-Perturbative Results 

In section 4, we considered the full non-abelian DBI action in the so-called adapted sym- 
metrized trace prescription. We showed how this can be motivated from the symmetrization 
procedure, and it can also be motivated as it reflects the regularity of the DBI action and 
the geometry of the probe brane embedding. As such, it gives physically sensible results 
even though the details of the quantitative results may not be exact and there may be 
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some other artifacts in special cases. Also for the unbroken phase this procedure is exact, 
as there is only one non-trivial generator in the DBI action. 

Overall, we found in 4.4, that the T — > divergences of the a' 2 expansion, e.g. in the 
free energy, are cured from the full DBI action. Inspecting the entropy, it turns out that 
the residual T = entropy of the unbroken phase vanishes in any of the broken phases, 
i.e. the ground state degeneracy gets lifted. Another effect of the regularity of the action 
is also that the critical temperature of the smaller- T c condensates gets lowered further, 
such they the critical temperatures are approximately exponentially separated. However 
the ordering of the critical temperatures does not change. 

While the p-wave phase that we discussed in sec. 4.1 does not change in its qualitative 
shape, the s-wave condensates do. Firstly, in all but the largest- T c phases, the second order 
phase transition changes to a first order one, and one of the compact s-waves crosses with 
the n = p-wave phase, and potentially the same happens at the higher modes. Also, the 
s-wave phases looked similar to each other and showed some multivaluedness that is usual 
for Minkowski - black hole ("meson melting") phase transitions. In the "flat" scalar in 
section 4.3, this was however a surprise. Still, we saw there the same, albeit now heavily 
distorted, attractor behavior of the internal parameters as in the a' 2 expansion. 

For the compact scalar, we found that the black hole embeddings are dis-favored com- 
pared to the unbroken phase for all but the largest T c cases. Their embeddings reach 
eventually a critical point at which the flux through the horizon to the probe brane (the 
un-condensated density) vanishes, and after which they do not extend to the horizon. In 
those "Minkowski" embeddings, the isospin flux/density is almost entirely generated by the 
condensate (the D5-D5 strings) in addition to some spread-out D3-D5 strings - in contrast 
to the abelian flux, where such embeddings would not be allowed [52, 41]. We found a 
remnant of the usual phase transition between those phases, but physically only a first 
order phase transition from the unbroken phase to the Minkowski phase takes place. 

For the largest-T c phase there exists still the second order phase transition to the 
black hole phase. But rather than seeing the expected first order phase transition to the 
Minkowski phase that would extend to T = 0, the gap in the free energy of the Minkowski 
phase and the value of the condensate diverge at a finite temperature. As shown in app. 
A. 2, at this point also the chemical potential vanishes and in the grand canonical ensemble 
there would be a broken phase in the UV limit. This unusual divergence is perfectly well 
understood from the gravity side as we discuss in section 4.2.3 and it is purely due to 
the fact that only a finite isospin flux is required to support these embeddings. Hence it 
is likely to be an artifact for the following reason: If we consider the symmetrization of a 
term that contains only one generator, e.g. ririnri, the adapted procedure is exact for this 
term. Hence, the overall "saturation" effects will persist if one takes the full DBI action. 
For very much mixed terms, like T1T2T1T2, however there will be some cancellation due to 
the symmetrization, e.g. of t\TiT\t<i and T2T1T1T2, and the result of the symmetrized trace 
will be lower than the "adapted" version. This implies that we over-estimate the coupling 
of the symmetry breaking field to the isospin potential in particular in the region where 
the fields are of the same order. This applies in particular to this case where we find the 
divergence. As this happens only when there is a particular balance of the coupling term, 
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we can expect that with a smaller coupling, the diverging phase would extend to large 
dimensionless density, i.e. vanishing temperature, where the divergence gets canceled by 
the overall density and temperature scaling, just as in the other phases. 

5.3 Outlook 

For future research, there may be three interesting avenues to take: Firstly, one could 
obviously consider more complicated configurations with more simultaneous fields. One 
particularly interesting one is the case where we explicitly separate the probe branes by 
turning on the T3-valued scalar in the flat directions. This corresponds in the field theory 
to separating the different layers. Hence it is a very physical problem to address. It does 
not affect the unbroken phase, but in the broken phases, a quick look at the DBI action 
suggests that it will lower the critical temperatures. Most interestingly, it will affect the 
p-wave and s-wave condensates differently. In particular in the simplest setup it would 
couple to all the fields except for the scalar in the flat directions. 

Another interesting thing to do would be to try to do a resummation of the expansion 
of the square root of the non-abelian DBI action after taking the symmetrized trace. In 
principle, this is unlikely to be successful, but this SU (2) case may be special, as we can 
use the straightforward anticommutators of the Pauli matrices. 
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A. Grand Canonical Ensemble 

To give one more motivation for the choice of the canonical ensemble, for completeness and 
because it may be interesting in some contexts, we show some results related to the grand 
canonical ensemble. 

A.l Perturbative Case 

We show the Gibbs free energy as a function of the "grand canonical temperature" 1/fl, 
in fig. 25. There we see that the behavior of the C2 — 0, n = mode is different from 
the other modes. Further, the energy gap is negative, so no symmetry breaking phase 
transition takes place and we see also no attractor behavior at small temperatures. Also, 
the first derivative is not continuous such that if there were a phase transition it would be 
first order. Still, this seems to be unphysical as the charge density in a superconductor (as 
opposed to e.g. a semiconductor) should be fixed. 
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Figure 25: The gap in the Gibbs free energy for the different condensates. As the gap is negative, 
we show the absolute value in the double logarithmic plot (right) and because of the strong diver- 
gence and large critical temperature of the (2 = 0, n = mode, we do not show this mode on the 
linear plot (left). 

A. 2 Non-Perturbative Case 
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Figure 26: Gap in the chemical potential A/x = ^unbroken — ^broken as a function of density at 
fixed temperature. The narrow black curve shows the chemical potential of the unbroken phase; 
the other lines are the same as in the plot for the free energy. Left: Double-logarithmic plot. Right: 
Logarithmic-linear plot. The inset shows the details of the £0 s-wave phase. 

Let us first look at the gap in the chemical potential Ap = p un broken — ^broken as a func- 
tion of density at fixed temperature, in fig. 26. There are three interesting observations: 
Firstly, in the limit of large densities, we see that A/2 is proportional to y/p, such that the 
chemical potential is lowed by an approximately constant fraction. Secondly, the gap near 
the "turning point" of the s-wave phases decreases with decreasing critical temperature 
(increasing critical p), and in the phase of <£ z A/2 is negative, i.e. the chemical potential 
increases, over a small range of density. This may be a feature of the regularity of the 
DBI action and corresponds also to the observation that the second order phase transitions 
turned into first oder ones. 

Thirdly, we see that the chemical potential actually vanishes at the "endpoint" of the 
(f2 = 0, n = phase. This can be perfectly well understood from the geometry, as the 
charge density p = — lim M _s.o d u p(u) remains finite, but the maximum uq becomes small (the 
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minimum radius of the probe brane becomes large) and hence p = hm u _>o p(u) = j u d u p(u) 
vanishes. 



Figure 27 : The gap in the Gibbs free energy 
of the different condensates as a function of the 
"grand canonical temperature" 1/p. The lines 
are the same as in the plot for the Helmholtz 
free energy, and the inset shows the details of 
o. 0.2 0.4 0.6 o.8 l. i.2 i.4 the small-temperature phases. 

This is reflected in the plot of the 27 where we show the Gibbs free energy as a function 
of 1/p, that serves as a temperature at fixed chemical potential. There, we see that the 
if 2 = 0, n = phase extends out to infinite temperatures and 5Q = ^unbroken — ^broken even 
becomes positive above some critical temperature. This feature seems to be an unphysical 
artifact of the adapted symmetrized trace prescription and should disappear if one takes 
the rigorous symmetrized trace and possible higher order corrections, as we will discuss in 
the conclusions. 

Else than that, the energy gap is negative in all the phases, so none of them is 
realized at fixed chemical potential and we see the lowering of the "critical temperatures" 
compared to the a' 2 expansion. Further, all the s-wave phases change in a similar fashion 
to a shape similar to the plots of the condensates - in the case of the compact scalar with 
a point at which the Minkowski and black hole phases intersect, similar to the meson 
melting phase transition [47, 41] in the abelian case. 
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